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FOREWORD 


The  present  compilation  of  abstracts  contains  recent 
information  on  seismometrv,  computer  seismology,  explosion  seismology, 
and  general  information  on  seismology,  including  some  material  on  the 
Rumanian  seismological  stations  and  several  articles  on  seismological 
research  in  Communist  China. 
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1 . Seismometry 


Aranovich,  Z.  I.,  and  O.  Ye.  Starovoyt. 

Status  of  observationswith  the  SD-I  long- 
period  seismograph  in  the  USSR.  Fizika 
Zemli,  no.  1,  1975,  88-94. 

The  introduction  of  SD-1  seismographs  into  base  stations 
of  the  Unified  System  of  Seismic  Observations  was  initiated  in  19t>8  as  part 
of  studies  of  the  upper  mantle  in  Europe  from  observations  of  long-period 
seismographs.  The  planned  network  of  stations  equipped  with  SD-1  seismo- 
graphs is  shown  in  Fig.  1.  At  the  present  time,  15  SD-1  seismographs 
operate  continuously  (see  Table  1) 


Fig.  1.  Location  of  seismic  stations  equipped 
with  SD-1  seismographs. 

1-  >perating;  2-  nonoperating  or  irregularly 
operating. 


Table  1 


Seismic  seismographs  stations  equipped  with  SD- 1 
(through  the  first  half  of  1972) 


No. 

Station 

Beginning  of 
operation 

Magnification 
N-S  E-W  Z 

Note 

1. 

Apatity 

Aug.  1 ' » 70 

- 

- 

970 

2. 

Bakuriani 

Dec.  1971 

- 

- 

560 

From  June  1972 
Vz  1400. 

3. 

Garm 

Jan.  '9hH 

100 

100 

- 

bi  1968-70 
magnification  was 
800. 

4. 

Irkutsk 

March  1971 

200 

200 

200 

5. 

Iul'tiu 

June  1970 

1000 

080 

980 

6 . 

Kishinev 

Jan.  1970 

1080 

In  Oct.  1972 
horizontal  compo- 
nents with  magni- 
fication 750  were 
installed. 

7. 

L' vov 

July  1969 

- 

- 

900 

In  July  1972 
horizontal  compo- 
nents installed 

8. 

No.  3 

Jan.  1971 

770 

770 

750 

Jam  1971  - Jan.  1972 
only  horizontal 
components . 

9. 

Moskva 

Jan.  1968 

- 

- 

1070 

10. 

Obninsk  SD-2(1) 
Press  - Ewing 

Jan.  196  7 
July  1966 

830 

800 

830 

800 

830 

1900 

11. 

Petropavlovsk 

May  1969 

1 100 

1100 

1050 

F rom  July  1972 
magnifications 
are  860,  960,  960. 

12. 

Pulkovo 

Feb.  1970 

- 

- 

100 

13. 

Simfe  rnpol ' 

Nov.  19  70 

830 

830 

980 

From  Nov.  1969 
only  vertical 
component. 

14. 

Sochi 

Aug.  D*71 

- 

- 

(1  500) 

15. 

Yuzhno  - 
Sakhalinsk 

Feb.  196  7 

1000 

1000 

1000 

Press  - Ewing 
seismograph 

-2- 


The  SD-1  seismograph,  consisting  of  SKD  seismometer 
(T  = 25  sec)  and  SPG  galvanometer  (T  f = 80-100  6ec),  is  characterized 
by  a passband  of  15-60  sec  at  the  0.  9 level  and  maximum  magnification  of 

1500.  It 6 threshold  sensitivity  at  V = 700-1000  was  determined  to  be 

SD  1 m o 

= 6 for  epicentral  distances  ^ = 20-100  . Its  dynamic  range  is 
^vmin  o 

not  sufficient  for  recording  PV  and  mpy  = 7 (A  = 20-40  ).  At  least  three 

sources  of  long-period  noise  were  identified  during  operation  of  SD-1 

seismographs  at  different  stations:  wind,  industrial  and  civilian  activity, 

and  fluctuations  in  atmospheric  pressure  and  temperature. 

The  records  of  the  same  earthquakes  made  by  SD-1  seismo- 
graphs installed  at  different  stations  are  characterized  by. 

1.  High  stability  of  waveform  of  body  waves  for  stations  with 
similar  epicentral  distance  and  azimuth  (Fig.  2); 

2.  Stable  pattern  of  dispersion  of  R waves  within  a period  range 
from  60  to  15-20  sac  but  with  amplitude  modulation  (Fig.  3).  Fig.  4 illustrates 
the  difference  in  amplitude  of  L surface  waves  from  the  earthquake  of 

5 May  1972  originating  in  the  Solomon  Islands  recorded  at  different  stations, 
which  is  suggested  to  be  due  to  high  directivity  of  seismic  radiation. 

3.  Distinguishable  R^  surface  waves  from  earthquakes  with 
M = 6.  5 - 6.  75  (Fig.  5). 

4.  High  stability  of  rripy  evaluated  frcm  these  records  (Table  2). 

The  fully  operational  SD-1  seismographs  provide  for  study  of  the  earth 
structure  along  the  following  profiles  and  triangles: 

1.  a.  Apatity -Pulkovo -No.  3- Kishinev  (Aleutions,  Balkans); 

b.  Apatity -Pulkovo -Obninsk -Simferopol'  (Alaska,  Turkey,  Red  Sea, 

East  Africa); 

c.  Apatity -Obninsk-Sochi  (Vancouver,  Tuamotu  Archipelago,  Red  Sea, 

East  Turkey,  Madagascar); 

d.  Kishinev-Simferopol'  -Sochi -Bakuriani  (Azores,  Central  America,  Iran, 

West  Indonesia,  New  Zealand). 


mr 
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Fig.  4.  Records  of  L surface  waves  from  the 
earthquake  of  5 May  1972  in  the  Solomon  Islands 
region,  made  by  SD-1  seismographs. 


Obninsk 

^ ^ Simferopol1 


Fig.  5.  Records  of  R^  surface  waves  from  the 
earthquakes  of  10  April  1972  in  Iran,  made  by 
SD-1  seismographs. 


Table  2 

Earthquake  magnitudes  determined  from  SD-1  records 


PH 

6.6 


o 

Earthquake 

Static  n 

A* 

mpv 

10  April  1972,  Iran 

Garm 

18 

- 

(0  = 28.  5°N, 

Sochi 

19 

6.  3 

> = 52.  8°E, 

Simferopol' 

22 

6 . 6 

mpV  = 6.  8, 

Kishinev 

26 

6.  5 

M = 6. 8) 

Obninsk 

29 

6.  7 

No.  3 

31 

6 . 6 

Pulkovo 

35 

6 . 6 

Apatity 

40 

6 . 6 

Irkutsk 

45 

6.  7 

Yuzhno -Sa - 
khalinsk 

69 

6.  6 

Iul'tin 

76 

6.6 

Petropavlovsk 

™ . 

6. 6 

6. 

6+0. 05 

1 1 April  1972, 

Kishinev 

67 

6. 4 

North  Atlantic 

No.  3 

69 

6 . 6 

ridge  (<[*  =1. 2°H, 

Simferopol' 

70 

6.  5 

X = 28.  4°W, 

Pulkovo 

74 

6.  5 

nipv  = 6.8, 

Sochi 

74. 

6.6 

M = 6.6) 

Obninsk 

75 

6 . 6 

Apatity 

78 

6.  4 

6.  5+0. 05 


TPV> 

tph* 

Note 

sec 

sec 

- 

14 

From  data 

22 

” 

of  5 stations 

17 

15 

Tot.  =26+4sec 

SH  - 

18 

16 

16 

^SH^PV^  1,63  | 

16 

- 

' 

16 

- 

i 

16 

- 

It 

- 

I 

12 

- 

14 

10 

14 

15 

16+2 

14+2 

14 

. 

; 

14 

- 

14 

- 

; 

13 

- 

i 

12 

- 

! 

12 

- 

1 

13 

- 

13+1 

' 

' 
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2.  a.  No.  3 -Obninsk -Pulkovo  (600j_50); 

b.  L'vov-No.  3-Kishinev  (6CKH1  50); 

c.  Kishinev-No.  3-Obninsk  (800  + 200); 

d.  Simferopol' -Obninsk-Kishinev  (850+350); 

e.  Simferopol1 -Obninsk-L'vov  (1000_+100). 


Feofilaktov,  V.  D.  , and  V.  A.  Masloboyev. 
Comparison  of  dynamic  characteristics  of 
standard  high-sensiti vity  seismographs  used 
in  the  USA  and  USSR.  Fizika  Zemli,  no.  1, 
1975,  95-100. 


A comparative  analysis  was  made  of  systematic  distortion  of 
the  shape  of  the  records  of  P waves  introduced  by  narrow-band  SKM-I1I  and 
Benioff  seismographs.  The  effect  of  self-inductance  on  the  motion  of  the 
narrow-band  Benioff  seismograph  was  considered  and  an  analytical  expression 
for  its  response  to  simple  harmonic  earth  motion  was  developed  by  the  method 
of  electromechanical  analogy. 


The  analog  circuit  for  Benioff  seismograph  used  in  the  analysis 

2 

is  shown  in  Fig.  1.  An  analysis  ci  the  formula  shows  that  for  T-*0  U-*T  /T.  , 

3 

and  for  T-*oo  U — ■ 1/sT  . The  response  curve  of  the  Benioff  seismograph 
calculated  bv  the  derived  formula  is  shown  in  Fig.  2. 


The  comparative  analysis  was  performed  in  terms  of  A = lg  A/T 
corresponding  to  maximum  pulses,  using  records  of  196  earthquakes  (^  = 
11-103°;  h from  normal  tj  530  km;  Mpy  = 4.  47-7.  3)  made  at  the  Obninsk 
station  in  1970.  The  apparent  periods  ranged  from  0.  4-2.0  sec. 


Fig.  1.  Electric  analog  for  Benioff 
seismograph. 


i-  analog  to  inertial  force  jU)m 

E - angular  velocity  0 of  the  galvanometer 

^ frame;  _ 

R = R + R.,  R = R + R,; 
s s 1 g g 2 

L-  self-inductance  of  the  seismometer  coil; 


1/m  - electromagnetic  converter  with  constant  m 

s D 

Gg-  magnetoelectric  converter  with  constant  Gg; 

Z,—  l/(/uini,  + 2i'.n"i . ♦ n , 


Zg  — 1 (/<iiA<  + 2f1,oAf  + A„//<u), 

where  e gQ  (egQ),  ns  (ng),  and  kg  are  air 

damping  coefficient,  angular  frequency,  and 
moment  of  inertia,  respectively. 


Fig.  2.  Response  curses  of  SKM-UI  and 
Benioff  seismographs. 

a-  calculated,  b and  c-  experimental. 
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The  variation  of  A A = asKM~  AB  with  aPParent  Periods 


. a ~ 

T is  shown  in  Fig.  3.  The  results  of  statistical  analysis  of  & A and 


SKM 

At  = t 


SKM  B 


-T  for  six  separate  period  intervals  are  given  in  Tables  1 and  2. 


The  variations  of  the  average  values  of  A A and  AT  with  the  average  values 


of  T„t,  calculated  for  the  six  period  intervals  are  shown  in  Fig.  4. 
SKM 


& fi 


Fig.  3. 


It  is  concluded  that  discrepancies  between  lg  A/T  and  T 
determined  from  the  records  of  SKM -III  and  Benioff  seismographs  are 
highest  in  the  period  range  0.6-1.  0 sec.  Namely,  values  of  lg  A/T 
determined  from  the  records  of  the  Benioff  seismograph  are  lower  by 
0.  1 than  the  corresponding  values  determined  from  the  records  of  SKM -III 
seismograph. 


+0,008 

0,000.1.3 

0.0230 


+ 0.050 
0,00256 
0,0505 

0.0152 

0.0481 


III 

IV 

25 

21 

+0,036 

+0.0!  1 

0,00380 

0.0GI5 

0.0123 

0,0344 

f ,076 


O.OOI56 

0,0396 


+0,003 

0.C0051 

0,0226 


+0,004 

0.00054 

0.0234 
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Nersesov,  i.  L.  , Ye.  I.  Gal'perin,  L.  M. 

Vorovskiy,  and  R.  M.  Gal'perina.  Alma-Ata 
high -sensitivity  radiotelemetric  seismograph 
test  site.  IN:  Fizika  Zemli,  no.  7,  1974,  72-77. 

A brief  description  is  given  of  the  set  up  and  instrumentation 
of  the  seismograph  stations  of  the  Alma-Ata  test  site,  as  well  as  the  initial 
results  of  its  operation. 

The  location  and  set  up  of  the  seimograph  stations  are 
illustrated  in  Fig.  1 and  Table  1.  The  network  consists  of  underground  stations 
set  up  in  deep  boreholes  (I,  V,  and  VII  in  Fig.  1)  and  surface  stations  (U,  III, 
IV,  and  VI).  The  underground  stations  are  equipped  with  SBU-V  seismometers 
with  natural  periods  of  1.  0 sec  and  preamplifiers  v/ith  negative  feed-back 
having  flat  temperature  curves  up  to  45°  C.  The  surface  stations  are  equipped 
with  SM-2M  seismometers  with  natural  frequencies  of  1 Hz  and  the  srjne  type 
of  preampliflier s.  The  signals  are  transmitted  to  the  central  station  at 
Novo  - Aleksey evka  (VII  in  Fig.  1).  The  power  of  the  transmitters  (with  power 
amplifier  units)  is  25  w;  without,  it  is  2.  5 w.  The  transmission  of  information 
via  RRS-IM  radio  relay  stations  is  based  on  double  conversion  FM-FM. 

Signals  are  recorded  on  the  same  strip  chart  by  an  RVZ-T  visual  recording 
system.  They  are  recorded  on  magnetic  tape  as  well  as  coded  using  an 
Udar-IM  analog -to -digital  converter. 

The  reading  accuracy  for  time  differences  in  0.  2 sec  in  the 
case  of  visual  recording  and  + 0.  03-0.  05  for  magnetic  recording.  Epicentral 
coordinates  are  determined  using  computer-generated  nomograms  for  given 
focal  depths  (see  Fig.  2). 
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Fig.  2.  Sample  Nomograms  Used  for  Determining 
Epicentral  Coordinates. 


^ - fir st -arrival  time  at  Alma-Ata  station; 
at  Talgar; 
at  Ozero; 

epicentral  distance  from  the  Alma-Ata  station 
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Rykov,  A.  V.  Calculation  of  the  main  parameters 
of  a magnetoelectric  seismograph.  IN:  Fizika 
Zemli,  no.  7,  1974,  78-83. 

A method  was  described  for  calculation  of  the  parameters 
of  a seismometer-galvanometer  system  based  on  the  feedback  ratio  concept. 
The  calculated  parametes  were  used  for  determining  resistance  in  seismo- 
graph circuits. 


In  the  calculation  of  the  maximum  magnification  V of  a 

max 

seismograph  for  the  case  of  a given  response  curve,  the  curve  is  specified 
graphically  in  Fig.  1.  Accordingly,  the  coefficients  in  the  response  equation 

i i|(i7  1 + u i fir1  + ^ r*  ■ % ... 


are  found  from  the  system  of  equations 

IKTj)  = 0.7;  U(T2)  = 1.0;  U(T3)  = 1.0;  U(T4)  = 0.7;  U(T5  ) = 0.  15. 


Fig.  1.  Pre-se*  response  curves  for  an 
SKM-3  seismograph  of  the  standard  type  III 
(1)  and  IV  (2). 


The  maximum  magnification  is  then  calculated  using  formula 

24  i j A . 

I 01»«  ~ 1/  “ 

ni  r a,  * 


where  A is  the  optical  arm  of  the  galvanomet  er. 


1 is  the  reduced  pendulum  length  of  the  seismometer,  kg  and  are 
inertia  moments,  TJ  is  the  bandwidth  factor,  and  kmax  is  feedback  ratio. 

In  calculating  seismograph  parameters  in  the  case  of  a given 
magnification  level,  the  desired  magnification  is  expressed  in  terms  of  the 
magnification  factor  F = £ I 'max'  w^ere  £ = ^^max  *8  t*ie  reduction 

factor  and  Fmax  = k/Tg2. 

The  table  below  gives  the  parameters  of  an  SKM-3  seismo- 
graph, standard  type  IV,  for  various  magnifications  and  the  parameters  of 
an  SKD  seismograph  for  the  maximum  and  standard  magnification. 


Table  1 


v to* 


SKM-3 


7.05 

0,382 

2,208 

0,683 

0,314 

1,561 

o:o 

0.324 

1 .028 

0,711 

0,369 

1,812 

5.0 

0,283 

1 ,721 

0,701 

0,414 

2,014 

4.5 

0.25  2 

1 ,040 

0,086 

0,434 

2,108 

4.0 

0,215 

1 .579 

0.667 

0,45! 

2,188 

3,5 

0,  !74 

1 .331 

0,648 

0,464 

2,252 

3.0 

0.134 

1 .501 

0,631 

0,474 

2,302 

2.5 

0.0% 

1 ,477 

0,617 

0,482 

2,341 

2.0 

0,90" 

1 ,459 

0,605 

0,488 

2,372 

1 .5 

3.036 

1 .446 

0.595 

0,493 

2,394 

SKD 


1,0 

0,228 

28.3 

0.468 

1,038 

0.884 

0 1 

0,0031 

25.0 

0,500 

1,200 

8,000 

k 

T. 

Ti 
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Tokmakov,  V.  A.  OSP  seismometer  and  its 
experimental  exploitation.  Trudy  Instituta 
fiziki  Zemli  AN  SSSR,  no.  16,  1974,  147-152. 


The  OSP  seismometer,  developed  jointly  in  1968-69  by  the 
Institute  of  the  Physics  of  the  Earth  (Soviet  Academy  of  Sciences)  and  the 
Special  Design  Office  of  the  Institute  of  Geophysics  and  Engineering  Seismology 
(Armenian  Academy  of  Sciences)  is  described.  A general  view  of  the 
seismometer  is  shown  in  Fig.  1 and  its  cross  section,  in  Fig.  2.  The 
OSP  seismograph  circuits  are  shown  in  Fig.  3. 


Fig.  1.  General  View  of  the  Vertical  OSPV 
and  Horizontal  OSPG  Seismometers. 


1-  magnetic -circuit  casing; 

2-  permanent  magnet; 

3-  pole  piece;  4-  induction 
coil;  5-  springs;  6-  interior 
cover;  7-  main  cover; 

8-  4 -conductor  cable;  9-  base 
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Fig.  3.  Electric  Circuits  of  the  OSP 
Seismograph. 

- operating  coil;  Rgg-  calibrating  coil; 
Rg-  galvanometer;  R^-R^  - resistances. 


The  specifications  of  the  OSP  seismometer  are  as  follows: 


1.  Mass  of  suspended  system  M , kg 

s 

2.  Natural  frequency  f , H 

S 2 

3.  Resistance  of  operating  coil  R , Ohm 

4.  Sensitivity  of  operating  coil  TnBg>  V/m/sec 

5.  Sensitivity  of  calibrating  coil  M , V/m/sec 

8 8 

6.  Damping  D , in  fractions  of  critical  damping 

s 

7.  Limit  velocity,  m/sec 

8.  Water-tightness,  in  meters  of  water  column 

9.  Weight,  kg 


4-  10' 


1.  25 


High  damping  is  achieved  using  aluminium  wire  for  the  induction  coil. 


The  OSP  seismograph  is  used  mainly  for  recording  the 
velocity  and  acceleration  of  ground  motion.  For  the  former,  the  OR -TIT 
(5  Hz)  and  M002  (25  H^,  = 15)  galvanometers  are  used;  for  the  latter, 

GB-IV  (120  H^)  and  M002  (120,  400  Hz)  galvanometers  and  the  N-735  loop- 
type  [recorder]  are  used. 


Beskorovayny/,  V.  L. , V.  P.  Marukhnenko, 

G.  S.  Sokolov,  and  A.  M.  Yampol'skiv,  Gain 
control  for  selection  signala  for  a digital  seismic 
system.  Othrytiya,  izobreteniya,  promyshlennyye 
obraztsy,  tovarnyye  znaki,  no.  18,  1974,  129. 

(Author  Certificate  428328). 

A patent  has  been  issued  for  the  above,  containing  a cascade 
amplifier,  a logic  circuit  for  amplification  selection  for  each  selection 
signal,  and  switches  for  commutation  of  the  outputs  from  the  amplifier  stages. 
To  reduce  the  equipment,  a converter  of  selection  signal  polarity,  a coin- 
cidence circuit,  and  a trigger  with  voltage  output  are  series  connected 
between  the  outputs  of  the  amplifying  stages  and  the  logic  circuit  inputs. 


Beskorovaynyy,  V.  L,. , V.  P.  Marukhnenko, 

G.  S.  Sokolov,  and  A.  M.  Yampol'skiy.  Digital 
multichannel  seismic  system.  Otkrytiya, 
izobreteniya,  promyshlennyye  obraztsy, 
tovarnyye  znaki,  no.  18,  1974,  129.  (Author 
Certificate  428329). 

A patent  has  been  issued  for  the  above,  containing  preamps 
connected  through  a channel  commutator  (having  a control  circuit)  to  the  main 
amplifier  which  has  digital  gain  control.  In  order  to  increase  the  dynamic 
range,  each  channel  in  the  system  has  a series  connected  delay  line  and  an 
additional  amplifier  with  digital  gain  control  connected  between  the  preamp 
output  and  the  channel  commutator  input.  The  additional  amplifier  is 
connected  to  the  output  of  an  additional  control  circuit  whose  inputs  are 
connected  to  at  least  two  sections  of  the  delay  line  and  the  channel 
commutator  control  circuit. 
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Slutskovskiy,  A.  I. , S.  V.  Nikiforov,  A.  I. 

Vasik,  and  V.  I.  Korinryy.  Method  of 
recording  ground  motion.  Otkrytiya,  izobreteniya, 
promyshlennyye  obraztsy,  £ov.*rnyye  znaki, 
no.  18,  1974,  129.  (Author  Certificate  428330) 

A patent  has  been  issued  for  the  above  method  which  is 
based  on  the  exposure  of  photosensitive  material  to  a light  beam  from  a 
controllable  source.  To  improve  visual  correlation,  the  signal  being 
recorded  is  split  into  two  channels.  In  the  first  channel,  the  signal  is 
converted  into  a voltage  corresponding  to  one  type  of  recording,  while  the 
signal  in  the  second  channel  modulates  the  recording  brightness.  In  this 
method,  the  signal  shape  is  changed  in  the  second  channel,  for  example, 
it  is  optimized,  high-frequency  modulated,  or  clipped. 
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2.  Computer  Seismology 


Gerver,  M.  L.  , V.  I.  Keylis -Borok,  Yu.  A. 

Kolesnikov,  A.  L.  Levshin,  V.  M.  Markushevich, 

B.  M.  Naymark,  and  V.  F.  Pisarenko.  Problems 
of  global  computational  geophysics.  Fizika  Zemli, 
no.  10,  1974,  33-45. 

A review  is  given  of  the  work  in  global  computational 
geophysics  performed  by  researchers  of  the  Institute  of  Physics  of  the 
Earth,  the  Institute  of  Applied  Mathematics  and  Chemical  Physics  of  the 
USSR  Academy  of  Sciences,  and  the  Moscow  and  Leningrad  State  Universities 
since  I960.  The  problems  considered  included  automation  of  primary  analysis, 
the  inverse  seismological  problem,  mathematical  modeling,  recognition 
algorithms,  and  evaluation  of  seismic  hazard. 


The  results  achieved  are  illustrated  in  Figs.  1-7. 

l’%  11.06.7Z  , J 66°,  7 1319  hrs 

30\  • 111 


Fig.  1.  Spectrum  of  a strong  earthquake  recorded 
by  the  SPS-300  seismograph  at  the  Obninsk  station. 

Length  of  record  - 10-16  hours; 
number  of  readings  - 10-16  thousand. 

Sharp  peaks  are  identified  with  natural  frequencies 
of  spheroidal  oscillations  of  the  Earth. 
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Fig.  2.  Seismogram  and  spectrum 
of  S waves  from  the  destructive 
Kamchatka  earthquake  of  24  Nov. 
1971  recorded  by  digital  equipment 
at  the  Naryn  station. 

High  energy  is  concentrated  in  the 
50-100  sec  range. 


Fig.  3.  Surface  waves  from  the 
destructive  Kamchatka  earthquake 
of  24  Nov.  1971  studied  using  the 
frequency -time  method. 

Group  velocity  of  Rayleigh  waves 
is  reliably  determined. 
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Fig.  4.  Use  of  the  frequency -time  analysis  in 
study  of  hydromagnetic  waves  (3  December  1965, 
Kerguelen).  A hydromagnetic  wave  is  split 
into  two  waves  shifted  in  phase. 


Fig.  5.  Determination  of  velocity  distribution 
from  time-distance  curves  for  refracted  waves 
from  a surface  source 


a-  in  the  case  of  two  waveguides;  b-  u(y) 
limited  by  b(y);  c-  additional  time-distance 
curves  for  refracted  waves  from  two  deep- 
seated  sources  between  and  below  waveguides. 
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Fig.  6.  Gutenberg  curve  g(r)  and  one  of 
"concurrent”  solutions  k(r). 

Cross  hatching-  velocity  distribution  in  the 
waveguide  v(r)  > b(r). 


Fig.  7.  Velocity  distributions  inferred  from 
data  from  deep  seismic  sounding  in  the  Central 
Turkmen  to  1962-64. 

1 - no  waveguide  is  assumed;  2-  a waveguide  with 
velocity  to  6 km /sec;  3-  a waveguide  with  velocity 
to  S.  6 km/sec. 
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Lander,  A.  V.,  A.  L.  Levshin,  V.  F.  Pisarenko 
G.  A.  Pogrebinskiy,  and  O.  Ye.  Starovoyt. 

I i * ■ 

Extraction  of  natural  oscillations  of  the  Earth  from 
records  of  the  Obninsk  observatory.  IN: 
Vychislitel'naya  seysmologiya,  no.  7,  Moskva, 
Nauka,  1974,  315-331. 


The  data  on  earthquakes  used  for  calculations  of  natural 
oscillations  of  the  Earth  are  given  in  Table  1.  The  frequency  response  of 
the  long-period  vertical  seismograph  used  for  recording  is  shown  in  Fig.  1 


Fig.  1.  Magnification  curves  of  vertical 
seismograph  SDS-300  (1)  and  Press-Ewing  (2). 

The  calculated  values  of  periods  of  identified  fundamental 
S (12  < n < 50)  and  overtone  . S (k>0)  oscillations  of  the  Earth  are  given 
in  Tables  2 and  3,  respectively;  values  of  quality  factor  Q are  in  Fig.  2. 
Dispersion  curves  determined  from  data  on  travelling  waves  are  shown 

in  Fig.  3. 


/♦  JO  Jf  Ht  M n 

Fig.  2.  Qp  for  fundamental  spheroidal 
oscillations  from  standing  (1)  and  travelling 
(2)  waves,  and  according  to  Anderson,  1967  (3). 


Data  on  earthquakes 
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Fig.  3.  Dispersion  curves  determined  from  R_, 
R .(1),  R R,  (2)  and  according  to  Dziewonski, 
1771  (3).  b 


Rudnitskiy,  V.  P.  Spectral  structure  of 
seismic  oscillations  and  approximate 
methods  for  its  study.  Geofizicheskiy  sbomik, 
no.  59,  1974,  36-45. 

Spectral  analysis  of  seismic  oscillation  was  performed  and 
the  following  conclusions  reached: 

The  spectral  analysis  of  seismic  oscillations  generated  by  pulse 
sources  in  the  solid  earth  is  adequate  only  as  a tool  for  approximate  study 
(see  Fig.  1). 
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The  separation  of  higher  harmonics  is  efficiently  accomplished 
by  band -pass  filters.  Weighting  of  a filtered  function  by  the  unfiltered  one 
improves  the  estimate  of  the  square  of  the  filtered  function  and  can  be  used 
for  furthe  r tracing  of  current  values  of  the  average  energy  spectra  along 
seismic  profiles. 

If  quasiharmonic  oscillations  are  generated  in  the  source  and 
the  condition  2Aco^  < Wq  is  satisfied  (where  OJ  ^ i3  the  carrier  frequency  of 
the  source  signal,  2Aa;^  band  width),  an  approximate  identification  of  the 
pulse  response  of  an  inhomogeneous  medium  can  be  reduced  to  a series  of 
linear  problems  by  means  of  band-pass  filters. 


L 


Soldatov,  V.  N.  Some  problems  of  earthquake 
signal  extraction.  AN  I3elSSR.  Izvestiya,  no. 
6,  1974,  129-130. 


The  probability  of  correct  detection  of  seismic  signals  is 
studied.  The  study  reduces  to  an  analysis  of  the  signal -to-noise  ratio  at 
the  output  of  a receiver  in  the  form 
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M(co  ) and  D(C>/m)  are  the  first  and  second  moments  of  the  normal  distribution 


m 


of  p(<p  /m)-  a function  obtained  at  the  output  of  the  H filter  if  the  output  function 


f(t)  contains  signal  m(t);  a is  the  signal  amplitude;  En  is  dispersion  of  white 
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noise;  and 


p ““S3” 

K- V 

H A«oi  ’ 


a*  


Y*  = 
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,AwL 

Am* 


(2) 


The  function  p (a,  /3,  a,  4>,  y,  E ) was  computed  accurate  to 

-7  n 

10  . Working  characteristic  curves  were  plotted  which  facilitate  determining 

the  optimum  width  and  frequency  response  position  for  a band-pass  filter 

with  a specified  probability  of  correct  signal  detection. 


Molchan,  G.  M.  , and  V.  V.  Ratushnyy.  Methods 
of  filtering  long-period  signals  over  short  time 
intervals.  IN:  Vychislitel'naya  seysmologiya, 
no.  7,  Moskva,  Nauka,  1974,  130-160. 

A comprehensive  account  is  given  of  filtering  methods  using 
exponential  smoothing  filters  with  transfer  functions  K^(p)  = 1 + pT)  ^ which 
were  developed  for  marine  gravimetry.  The  accuracy  of  measuring  a smoothed 

function  at  random  points  is  analyzed.  Numerical  examples  are  given 

n •: 

illustrating  the  design  of  filters  in  the  form  K ^ T.  < Tq  and 

measurement  of  a smoothed  function  at  extrema  of  a centered  function  (l-K)s(L). 


It  is  concluded  that  efficient  low-frequency  filters  can  be 
designed,  combining  filters  of  exponential  smoothing  by  multiplication  and 
addition.  Measurement  of  the  function  at  zeros  of  an  auxiliary  function  with 
zero -average  value  can  be  efficient  only  if  the  measured  and  auxiliary 
functions  are  sufficiently  well  correlated.  The  highest  accuracy  is  attained 
if  zeros  of  the  auxiliary  function  are  determined  by  a vertical  window-type 
operator.  Measurements  at  random  points  over  long  time  intervals  are  equal 
to  measurements  at  uniformly  spaced  points  with  the  same  density  in  the  sense 
of  an  asymptotic  equality  of  depressions  of  the  average  error  of  all  measure- 
ments. 
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Borkovskiy,  G.  M.  , A.  A.  Drozdov,  V.  M. 
Kreysberg,  and  A.  A.  Tsukanov.  Device  for 
reproduction  of  digital  seismic  information.  Otkr, 


izobi,  no.  20,  1974,  115.  (Translation) 

The  device  described  by  this  patent  contains  a analog -code 
converter,  a dynamic -range  reproduction  circuit,  channel  selector, 
continuous  analog  signal  shaper,  and  recorder.  To  improve  economy  and 
provide  identical  channel  control,  between  the  output  of  the  dynamic  range 
reproduction  circuit  and  the  imput  of  the  channel  selector,  an  automatic 
gain  control  system  (AGC)  has  been  added  for  multichannel  seismic  information 
in  a single  amplifying  channel.  In  the  AGC,  the  output  of  the  reproduction  circuit 
is  connected  to  the  amplification  input  of  a common  amplifier  with  a regulated 
gain  coefficient,  whose  output  is  connected  to  the  input  of  the  amplifier  of  the 
AGC  loop  of  the  common  amplifier.  This  amplifier  is  connected  through  a 
series  connected  phase  shifter,  rectifier,  and  filter  to  the  control  input  of 
the  common  amplifier  which  is  connected  to  the  input  of  the  channel  selector. 

The  outputs  of  any  two  opposing  end  channels  and  of  any  central  channel 
of  the  continuous  analog  signal  shapers  are  connected  to  the  input  of  a 
mixer.  The  mixer  output  is  connected  to  the  amplifier  input  of  AGC  loop. 

To  the  second  inputs  of  the  phase  shifter,  a switchable  splitter  is  connected, 
the  first  input  of  which  is  connected  to  the  output  of  a multivibrator,  while 
the  second  input  of  the  splitter  is  connected  to  the  output  of  a time  relay. 

The  time-relay  input  is  connected  to  the  circuit  providing  the  detonation 
time  mark. 

Borkovskiy,  G.  M.,  A.  A.  Drozdov,  and 
V.  M.  Kreysberg.  Method  of  visualization  of 
multichannel  digital  seismic  information.  Otkr 
izobr,  no.  20,  1974,  115.  (Translation) 


Digital  information,  compressed  in  time,  is  converted  into 
analog  form  with  subsequent  channel  distribution.  This  method  is  distinguished 
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by  the  amplification  of  compressed  analog  seismic  signals  employed  for 
the  purpose  o'  compressing  the  dynamic  range  in  the  case  of  identical 
channels  and  co  preserve  dynamic  distinctiveness  of  the  records.  The 
amplification  coefficient  is  changed  proportionally  to  the  summed  energy 
of  the  signals. 


Gur'yanov,  V.  M.  Breakdown  of  a seismic  wave 
field  into  travelling  waves.  IN:  Sb.  Primeneniye 
TsVM  i sredstv  vychisl.  tekhn.  v geol.  i geofiz. 

Saratov,  1974,  3-18.  (RZhMekh,  11/74,  #11V837) 

(Translation) 

The  theory  of  the  breakdown  of  summary  seismic  wave  field 
into  its  component  travelling  plane  waves  was  developed  on  the  basis  of  the 
solution  of  the  Lame  equations  for  an  nonhomogeneous  isotropic  elastic  theory 
with  ray  approximation.  This  approach  is  associated  with  a study  of  character- 
istics of  spectral  functions  which  are  obtained  by  applying  a two-sided  integral 
Fourier  transform  to  the  expression  for  the  displacement  field  for  an  non- 
homogeneous elastic  medium.  It  was  shown  that  the  amplitude -frequency 
characteristics  of  seismic  waves  recorded  by  a linear  receiver  arrangement 
satisfy  a system  of  integral  equations  of  the  Fredholm  type,  which  can  be  solved 
by  an  approximate  method.  Coefficients  of  the  approximate  solutions  are 
found  from  the  condition  of  functional  minimization  on  the  basis  of  the  least 
squares  method. 


The  basic  principles  of  the  theory  developed  are  illustrated  by 
an  example  of  an  approximate  breakdown  of  summary  seismic  field  into  waves 
with  equal  time-distance  curves.  It  was  pointed  out  that  it  is  sufficient  to 
have  a minimum  volume  of  experimental  information:  time -distance  curve 
for  first  arrivals  and  time  distance  curve  for  trailing  edges;  information  on 
the  shape  and  intensity  of  interference  waves  is  not  needed. 
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Lander,  A.  V.  Methods  for  interpretation 
of  results  of  frequency -time  analysis.  IN: 

Vychislitel'naya  seysmologiya,  no.  7,  Moskva, 

Nauka,  1974,  279-315. 

A comprehensive  account  of  frequency -time  analysis  of 
seismological  data,  fundamentals  of  the  analysis  of  body  and  surface  seismic 
waves,  as  well  as  a description  of  the  analysis  techniques  are  given.  The 
following  is  a summary  of  the  study. 

The  physical  meaning  of  the  results  of  frequency -time  analysis 
(FT A)  depends  strongly  both  on  the  shape  of  the  signal  and  the  techniques  of  the 
analysis  applied.  It  varies  widely  depending  on  the  selection  of  filters. 
Therefore,  the  selection  of  filters  is  the  most  important  step  in  FTA.  The 
main  problems  - extraction  of  signal  and  measurement  of  its  parameters  are 
solved  simultaneously  by  FTA.  It  is  not  an  optimum  method,  for  the 
sign j.1  extraction  of  parameter  measurements,  and  it  has  limited  possibilities. 
The  best  signal  extraction  may  be  accompanied  by  poor  parameter  measurements 
and  vice  versa.  In  every  real  situation,  a compromise  should  be  selected.  The 
selection  of  filters  is  determined  by  the  problem  to  be  solved,  as  well  as  by 
the  type  of  signals  studied.  A possible  way  to  improve  methods  of  FTA  is  a 
preliminary  transformation  of  the  signal  into  a form  for  which  the  conflict 
between  the  extraction  of  a signal  and  measurement  of  its  parameters  is  less 
severe.  The  great  advantage  of  FTA  is  that  it  requires  no  prior  information 
on  the  signal. 
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3.  Explosion  Seismology 


Kostyuchenko,  V.  N. , and  V.  N.  Rodionov. 
Radiation  of  seismic  waves  during  powerful 
underground  explosions  in  solid  rock. 
Fizika  Zemli,  no.  10,  1974,  65-73. 


An  estimate  is  made  of  parameters  of  seismic  waves 
radiated  by  underground  explosions,  based  on  an  approximate  model 
of  underground  explosions  in  solid  rock  as  developed  by  Radionov  et 
al. , 1968  and  Radionov  et  al. , 1971.  The  estimates  are  compared  to 
experimental  data  on  underground  nuclear  explosions  in  rock  salt  and 
granite. 


The  results  of  calculations  of  the  maximum  displacement 
and  velocity  of  comp*-essional  waves  from  explosions  in  rock  salt  and 
granite  are  shown  in  Figs.  1 and  2. 


Fig.  1.  Maximum  displacement  (a)  and  velocity  (b)  in  compressional 
waves  from  explosions  in  salt. 

Curves  1,  2,  and  3-  calculations  for  absorption  decrement  D = 0,  0.  05, 
and  0.1,  respectively;  Curve  4-  average  from  measurements  in  the  interior 
and  at  the  s irface  of  the  rock  mass  (open  and  solid  circles,  respectively); 
Curve  5 and  crosses  - residual  displacements  from  Salmon  explosion, 
(circles  - Gnome  data;  triangles  - Salmon;  squares  - USSR) 
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Fig.  2.  Maximum  displacement  (a)  and  velocity  (b)  in  congressional 
waves  from  explosions  in  granite. 

Curves  1,  2,  3-  calculations  for  D = 0,  0.05,  and  0.1,  respectively: 
Curve  4-  average  from  measurements  in  the  interior  of  the  rock  mass 
(open  circles  - Sahara  data;  triangles  - Hardhat  and  Shoal). 

The  maximum  radius  of  the  elastic  zone  a.^,  maximum  displacement  at 
the  boundary  of  the  elastic  zone  uQ,  and  characteristic  radiation  time  T 
we  re  ^calculated  to  be:  a /q1 ' 3 = 1 m/kg1/3,  uQ/q1/3  = 0. 1 cm/kg1 /3,  ° 
T„/q  = 1 msec/kp  for  salt  (^=  2.2  gr/cm3,  Vp  = 4.4  km/sec. 


*VV  P 


Tg/q  ^ 7 J - 1 m8ec/kgJ 

10  );  and  0.  87  m/kg1 /3,  8.  7 x 10'2  cm/kg1/3,  and  0.  64  msec/kg1/3 

for  granite  (p  = 2.  7 gr/cm3,  Vp  = 5.  5 km/sec,  a,/pVp2  = 10'3). 


Measured  and  calculated  time  dependencies  of  displacements 
from  explosions  in  salt  and  granite  are  shown  in  Fig.  3 and  Fig.  4,  a. 
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Fig.  3.  Measured  displacement  fields  from 
Salmon  explosion  (solid  lines)  and  calculated 
values  (dashed  lines). 


. kb»r 


Fig.  4.  Displacement  (a)  and  radial  stress  (b)  from  Hardhat  explosion  and 
calculated  values. 

a.  solid  line  - measurements;  dashed  line  - calculations. 

b.  thin  solid  line  - measurements  at  R = 120  m;  dashed  line  - calculations 
for  R = am  = 150  m according  to  gp  = 1/Rn  (n  = 1.  5-2);  thick  solid  line  - 
calculations  according  to  -gp  (r)  = <7*  (1-e-dr),  where  a = 2.2/tj^, 

g*  = 10-3pVp2  = 800  b. 
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it  i.s  pointed  out  that  a quasistatic  approximation  adopted 
in  the  analysis  of  the  last  stage  of  explosion  development  can  describe 
adequately  only  radiation  in  the  low-frequency  range,  as  illustrated  in 
Figs.  4 b and  5.  Figure  5 gives  the  spectrum  of 


Fig.  5.  Spectrum  of  signals  radiated  during 
explosions. 

Open  circles  - Gnome  data;  solid  circles  - 
Hardhat;  Curve  1 - calculations  in  quasistatic 
approximation;  Curve  2 - calculations  for  head 
waves;  Curve  3 - calculations  by  |^[m1qq  q‘  flm 

p or 

the  derivative  of  potential  |$>'l  which  at  large  distances  represents  the 

spectrum  of  displacement  in  compressional  waves.  Calculations  of  |<$>| 

were  made  assuming  y = 0.  5,  a a^/Vp  = ^ an<^  ’ )•  t,„. 

where  a = 2.  2/t  in  a quasistatic  approximation  and  cr  „ = 0,e  alr,  where 
m i w 

a - 10  Qt , a I = a „ for  head  waves.  The  figure  illustrates  that  for  oj  < 1-2, 
the  quasistatic  approximation  sufficiently  well  describes  the  spectral 
amplitude  of  waves,  while  at  high  frequencies,  calculations  for  head  waves 
give  better  results. 
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Zel  manov,  I.  L. , V.  N.  Kologrivov,  A.  A. 

Krasavin,  V.  1.  Kulikov,  V.  V.  Pedanov,  and 
A.  M.  Tikhomirov.  Study  of  seismic  effect 
of  an  explosion  on  a transparent  model.  Fizika 
Zemli,  no.  10,  1974,  80-91. 

A description  is  given  of  the  techniques,  method,  and  results 
of  a model  study  of  seismic  effects  of  explosions.  The  results  were 

compared  to  those  obtained  from  data  on  nuclear  and  chemical  explosions 
in  solid  rock. 


The  models  used  were  fabricated  from  K-8  and  TF-5  optical 
glasses  and  aluminum -potassium  alum  (KA^SO^  • 12HzO).  Mechanical 

and  physical  properties  of  materials  used  for  the  models  are  given  in 
Table  1. 


Table  1 


Material 

p, 

gr/cm3 

E’  7 
kg/ cm4 

V 

kmf  sec 

km/sec 

kg/cm2 

n 

K-8  glass 

2.  52 

DC 

►— « 
O 
en 

0.  19 

6.  06 

3.  74 

7-10. 103 

1.  516 

TF -5  glass 

4.  77 

5.  8- 105 

0.  22 

3.  71 

2.  22 

6.  5*  103 

1.  755 

Aluminum 

1.  75 

2.  25- 105 

0.  25 

3.93 

2.  26 

150 

1.456 

Not  e : Q 

is  compressive  strength. 

n refracti 

ve  index. 

The  explosions  were  simulated  by  a plasma  focus  produced 
by  focusing  single-pulse  radiaLon  of  a Q-switched  ruby  laser  in  the  models, 
as  shown  in  Fig.  1.  The  parameters  of  the  laser  are  as  follows:  pumping 
energy  3kj;  pulse  energy  lj;  pulse  half-width  20  nsec;  and  beam  diameter 
about  10  mm.  Energy  density  of  the  "laser"  explosion  was  .l.Sx  104  j/cm3. 
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Fig.  1.  Experimental  Set-Up. 

1-  laser;  2-  sloping  plate;  3-  power-level 
indicator;  4-  telescope;  5-  target;  6-  Mach- 
Zehnder  interferometer;  7-  liquid-dye  laser; 

8-  spherical  mirror;  9-  high-speed  camera. 

The  parameters  of  compres sional  waves  were  determined  by  the  interference 
method.  Interferograms  v/ere  processed  by  the  "travelling  wave"  method. 

An  example  of  interferogram  processing  of  the  explosions  in  TF-5  is  shown 

in  Fig.  2. 
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Fig.  2.  Displacement  of  the  central  band  (a)  density 
diagrams  (b)  and  second  derivatives  of  displacement 
potential  (c). 

Curves  1,  2,  and  3 represent  2.  27,  5.  95,  and  11.  60  mm 
offsets  of  the  camera  slit  from  the  explosion  center, 
respectively. 
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The  results  of  model  studies  of  explosions  in  solids,  as  well 
as  their  comparison  with  results  obtained  from  data  on  nuclear  and  chemical 
explosions  in  rock,  are  given  in  Figs.  3-6  and  Tables  2-4. 
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Fig.  3.  Displacement  potential  and  its 
derivatives  for  explosion  in  TF-5  (energy  1 J). 
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Fig.  4.  Displacement  potential  and  its 
derivatives  for  explosion  in  alum  (energy  1 J). 
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Fig.  5.  Diagrams  of  velocity 
and  displacement  in  TF-5  (a) 
and  alum  (b)  at  6 mm  from  the 
explosion  center. 


Fig.  6.  Maximum  mass  velocity 
for  laser  and  nuclear  explosions 
(1-  after  Derlich,  1970;  2-  Trembly 
and  Berg,  1966). 


Table  2 

Radius  of  shattered  and  fractured  zones 


Medium 

Explosion 

kg/cm^ 

rl 

m/kT1/3 

r2  1 
m /kT 

[ 

K-8  glass 

— — — - - - - 

Laser 

o 

2 

i 

r 

8.8 

20.  5 

TF-5  glass 

-11- 

6 500 

5.  74 

15.  8 

Alum 

-11- 

160 

39 

130 

Granite 

Nuclear 
(Derlich,  1970) 

2000 

10.2 

27 

(?) 

-11- 

(Kedrovskiy  et  al.  . 
1970) 

20-30 

| 

50-70 

Rock  salt 

Chemical 
(Rodionov  et  al. , 
1967) 

400-600 

33 

. j 

Sodim 

160-200 

44 

i 

thiosulphate 

-11- 

- 1 

.1 

1 


Seismic  yield  of  laser  and  nuclear  explosions 


Material 

Portion  of  explosion  yield  e,  °/c 

K-8 

] 

0.  09 

TF-5 

0.  12  1 

Alum 

0.92 

Tuff  (Rodionov  et  al.  , 1971) 

0.  2-0.  3 

Granite  (Trembly  and  Berg,  1966) 

0.  7-2. 0 

Table  4 


Ejecta  volume  for  laser  and  cavity  volume 
for  nuclear  and  chemical  explosions 


Material 

Explosion 

V,  m3/kT 

V,  m3/kT 

K-8  glass 

Laser 

28.  3 

- 

TF-5  glass 

I! 

38.  3 

Alum 

11 

1250 

Granite 

Nuclear  (Derlich, 
1970) 

- 

1640 

Granite 

" (Kedrovskiy 
et  al.  , 1970) 

- 

4000-14000 

Rock  salt 

" (Boardman, 
1970) 

- 

3500-7900 

Rock  salt 

" (Kedrovskiy 
et  al.  , 1970 

5600-100000 

Rock  salt 

Chemical  (Rodionov 
et  al. , 1967) 

- 

5450 

Sodium  thiosulphate 

If 

12700 

, - - j 

Vovk,  A.  A. , A.  V.  Mikhalyuk,  and  I.  V. 

Belinskiy.  Producing  underground  cavities 
by  a confined  explosion  in  an  isotropic  rock 
mass.  PM,  v.  10,  no.  9,  1974,  53-59. 

Formulas  are  developed  for  cavity  radius  produced  by  an  under- 
ground explosion,  taking  into  account  the  behavior  of  rocks  under  blast  load, 
which  are  suitable  for  underground  engineering  applications.  The  results  of 
experiments  with  uniaxial  dynamic  compression  of  rocks  without  lateral 
expansion  are  given  in  Fig.  1.  The  calculating  model  is  shown  in  Fig.  2. 


Fig.  1. 


The  results  of  calculation  and  their  comparison  with 
experimental  data  are  given  in  the  table. 
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Koryavov,  V.  P.  Forerunner  effect  on  motion 
in  the  near -explosion  zone.  Fizika  Zemli,  no. 

10,  1974,  74-79. 

The  results  are  presented  of  an  analysis  of  the  effect  of 
an  elastic  medium's  strength  on  elastic  energy,  ca\ity  dimension,  and 
fracture  zone  dimension,  resulting  from  a powerful  underground  explosion. 

In  the  formulation  of  the  problem,  elastic  waves  are  radiated  at  the  front 
of  destruction  of  an  elastic  medium,  i.  e. , its  conversion  into  a fluid  medium. 
The  fluid  medium  is  described  by  equations  of  state  obtained  from  shock 
adiabatics  for  solids  (Koryavov  and  Vilenskaya,  1968).  In  the  present 
calculations,  a change  in  the  conditions  at  the  shock  wave  front  is  introduced. 

All  calculations  are  made  for  fracture  condition  in  the  form 

pe  - Pe  =a- 

where  p and  p_  are  radial  and  azimuthal  stress,  respectively  in  the  elastic 
6 W 

wave  in  front  of  the  shock  wave.  It  is  assumed  that  a = 5 kb,  c = 5.  44  km /sec, 

j£ 

a = 0.  27,  and  other  parameters  are  the  same  as  in  the  above-mentioned  work. 

The  results  of  calculations  are  shown  in  Figs.  1-3.  A 
comparison  of  the  present  results  with  experimental  and  theoretical  values 
reported  by  other  authors  is  given  in  the  table. 


Fig.  1.  Time  variation  of  elastic  energy  outflow  at  shock 
front  (a)  and  propagation  of  shock  front,  forerunner,  and 
cavity  boundary  (b). 

Dashed  lines  - propagation  of  shock  front  and  cavity  boundary 
without  accounting  for  elastic  waves. 
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Fig.  2.  Mass  velocity  profiles  (D  = 3.7  km/sec) 
considering  (solid  lines)  and  not  considering  (dashed 
lines)  elastic  forerunner. 


Fig.  3.  Elastic  potential  and  its  derivatives 
reduced  to  a detonation  energy  of  1 kt. 
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4.  General 


Sinitsyn,  A.  P.  Seismic  effect  from  several 
sources.  AN  SSSR.  Trudy  Instituta  fiziki 
Zemli,  no.  16,  1974,  65-72. 

An  analysis  is  made  of  seismic  effect  due  to  a point  force 
suddenly  applied  to  the  surface  at  an  elastic  halfspace  in  the  case  of: 

1)  two  fixed  sources;  2)  two  moving  sources  approaching  each  other;  and 
3)  two  pairs  of  moving  sources,  two  approaching  each  other  and  two 
receding  from  each  other. 


The  wave  patterns  in  the  first  case  are  shown  in  Figs.  1 and  2. 
In  the  inner  zone  of  the  halfspace  (between  the  sources),  seismic  effect  is 
amplified  due  to  the  interaction  of  waves.  The  collision  area  (shaded  area 
in  Fig.  1)  can  be  considered  as  an  additional  source,  twice  as  strong  as  the 
individual  sources.  The  wave  pattern  in  the  outer  zone  is  shown  in  Fig.  2. 
The  seismic  effect  in  this  zone  is  amplified  due  to  the  "imaginary"  source 
in  the  inner  zone,  and  exceeds  the  simple  summary  effect  of  two  sources  by 
a factor  of  1.  5 (see  Fig.  3). 

In  the  second  case  which  is  illustrated  in  Fig.  4,  the  two 
sources  approach  each  other  at  supersonic  speed.  Due  to  interaction  of 
waves,  compressional  wave  reflection  occurs  which  approximates  reflection 
from  a fixed  obstacle.  Thus,  the  amplitude  of  reflected  waves  is  twice  as 
large  as  the  amplitude  of  incident  compressional  waves. 

In  the  third  case,  which  is  illustrated  in  Fig.  5,  the  wave 
Pa^ern  the  inner  zone  is  the  same  as  in  the  second  case.  In  the  outer 
zone,  three  different  compressional  waves  occur:  direct  waves  from  two 
pairs  of  sources  and  reflections  from  the  collison  plane  which  is  formed  in 
the  inner  zone.  The  seismic  effect  in  the  outer  zone  exceeds  the  simple 
summary  effect  of  two  pairs  of  sources  by  60-70T(  (see  Fig.  6) 
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Fig.  1.  Collision  of  waves 
from  two  sources. 


Fig.  2.  Generation  of  an 
imaginary  source. 


Fig.  3.  Diagram  of  wave  network 
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Fig.  5.  Waves  resulting  from  two 
pairs  of  moving  sources. 


Fig.  6.  Diagram  of  wave 
network. 


AN  SSSR.  Institut  fiziki  Zcmli.  Apparatura 
i metodika  seysmometricheskikh  nablyudeniy 
■ SSfc*f  ( Equipment  and  methods  for  scismomctric 
observations  in  the  USSR)  Moskva,  Izd-vo 
Nauka,  1974,  242  p. 


This  book  reviews  fundamentals  in  the  theory  of  seismometer 
channels,  and  the  basic  types  of  seismic,  recording,  and  control  instrumental 
used  in  permanent  observatories,  field  stations,  and  engineering  sites  are 
described.  Methods  and  equipment  used  at  permanent  and  temporary 
observatories  are  reviewed.  The  table  of  contents  of  this  book  is  as  follows: 
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Savarenskiy,  Ye.  F.  , and  G.  L.  Kosarev. 

Effect  of  crustal  structure  beneath  a 
seismological  station  on  motion  in  compres- 
sional  waves.  Fizika  Zemli,  no.  10,  1974, 

113-120. 

An  analysis  is  made  of  the  dependence  of  apparent  angles 
of  emergence  of  compressional  seismic  waves  at  a seismological  station 
on  local  crustal  structure,  real  angles  of  emergence  at  the  crustal  base 
and  frequency  of  oscillations  as  well.  An  attempt  is  made  to  determine 
the  parameters  of  a layered  crust  from  spectra  of  apparent  angles  of 
emergence. 


Computer  calculated  values  of  vertical  and  horizontal 
components  of  the  amplitude  of  the  motion  of  the  earth's  surface  as  recorded 
by  a long-period  seismograph  for  a single-layered  crust  overlaying  a half- 
space and  for  different  incident  pulses  are  shown  in  Figs.  1 and  2.  Observed 
and  calculated  seismograms  for  a multi-layered  crust  overlaying  a half- 
space are  shown  in  Fig.  3 and  4.  A comparison  between  seismograms  calculated 
for  a single-layered  crust  overlaying  a single-  and  multi-layered  upper 
mantle  is  given  in  Fig.  5.  The  model  of  a 1000  km-thick  multi-layered  upper 
mantle  (several  hundred  2-5  km-thick  layers)  used  in  the  calculations  was 
based  on  the  velocity  distribution  for  shear  waves  typical  of  continental 
platforms.  Densities  were  calculated  by  the  formula  pk  = 0.  763  + 0.  328  ak 
(where  a is  the  velocity  of  compressional  waves). 
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Fig.  2.  Motions  in  compressional  waves  at 
the  surface  of  a single-layered  crust  as  recorded 
by  a long-period  seismograph  for  different  periods 
of  incident  pulse. 
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Fig.  1.  Motions  in  compressional  waves  at  the 
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incident  pulse). 

a-  real  motions;  b-  record  by  a long-period 

seismograph  (Ts  = 30  sec;  T = 100  sec, 
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Fig.  3.  Calculated  (thin  line)  and  observed  (thick 
line)  seismograms  of  compressional  waves  from 
the  27  May,  1970  earthquake  in  the  Bonin  Islands 
and  observed  at  the  Obninsk  station. 


Fig.  4.  Calculated  (thin  line)  and  observed  (thick 
line)  seismograms  of  compressional  waves  from  the 
24  March,  1967  earthquake  in  the  Malay  archipelago 
and  observed  at  the  Yuzhno-Sakhalinsk  station. 
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the  surface  of  a single-layered  crust  overlaying 
a single-layered  (thin  line)  and  multi-layered 


Sobolev,  G.  A.  , and  L.  1).  Slavina.  Rapid 
changes  in  electric  and  seismic  properties 
of  the  medium  in  a seismically  active  region. 
DAN  SSSR,  v.  215,  no.  5,  1974,  1101-1104. 


Data  on  rapid  changes  in  electrotelluric  field  and  velocities 
of  compressional  and  seismic  shear  waves  from  earthquakes  withK  > 8, 

^ < 150  km  observed  in  the  Kamchatka  region  were  analyzed  using  uniform 
processing  techniques.  The  observation  stations  and  distribution  of 
epicenters  for  strong  earthquakes  are  shown  in  Fig.  I. 


Fig.  1.  Positions  of 
strong  earthquakes. 


bservation  stations  and  epicenters  of 
Shad  in  c indicates  that  prior  to  the  earth  - 

Station 
at 


quakes  there  were  observed  peaks  in:  1-  J£V|’  at  No.  1 St 
|5V|“  at  No.  3,  3-  |A  ' ' v | 2 at  No.  2;  4-  |A  v /v  | 2 

No.  3;  5.  2 at  NcV.  ^ P s 
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Variations  in  the  electrotelluric  and  elastic  terrestrial 
fields  in  the  Kamchatka  region  during  1972-73  are  illustrated  in  Fig.  2. 


Fig.  2.  Variations  of  at 

No.  1 and  No.  J stations  (1,  2) 
and  |av  /v  I 4 at  the  No.  2.  No. 
3,  and  Nb.  4 stations  (3,  4,  5). 


Earthquakes: 


I-  18  January  1972,  M = 5.4; 

II-  2 5 May  1972,  mpv  =6.2; 

III-  26  June  1972,  mpv  = 6.  1; 

IV-  2 August  1972,  M = 5.  9; 

V-  4 August  1972,  M =6.7; 

VI-  27  September  1972,  M = 5.  1; 

VII-  10  November  1972,  mpv  5.2; 

VIII-  25  December  1972,  M 5.6; 

IX-  25  January  1973,  M 5.0; 

X-  28  February  1973,  M - 7.4. 


It  was  found  that  numerous  strong  earthquakes  are  preceded  by 
rapid  changes  in  electrotelluric  field  and  the  ratio  between  velocities  of 
compressional  and  seismic  shear  waves.  I hcse  changes  were  not  strictly 
simultaneous  over  the  entire  region,  but  were  observed  at  rather  large 
distances  from  the  subsequent  earthquakes. 
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Yanovskaya,  T.  B.  Determining  depth  of  the 
core  boundary  from  travel  times  of  PcP  waves 
and  time -distance  curve  of  P waves.  IN:  Sb. 

Teor.  i vychisl.  geofiz.  No.  2.  M. , 1974,  76- 
90.  (RZhMekh,  11/74,  #11V828).  (Translation) 

A method  is  proposed  for  determining  the  radius  of  the 
earth's  core  and  the  velocity  gradient  in  the  layer  above  the  core,  using 
data  on  travel  times  of  PcP  waves  and  time-distance  curve  of  P waves. 

The  formulas  for  transformation  of  time-distance  curves  of  refracted 
waves  into  time -distance  curves  of  reflected  waves  were  used.  It  was 
shown  that  the  ambiguity  in  time -distance  curves  of  P waves  at  a distance 
of  about  20°  (loop,  shadow  zone)  has  practically  no  effect  on  time-distance 
curves  calculated  for  PcP  waves.  The  core  radius  r.  and  the  velocity  gradient 
dv/dz  were  estimated  using  previously  published  data  (see  Herrin  E.  , 

Tucker  W. , Taggart  J.,  Gordon  D.  W.  , Lobdell  J.  L. , Bull.  Seismol.  So*  . 
Amer. , 1968,  58,  No.  4,  1273-1291).  A 90 °Ji-  confidence  region  for 
calculated  parameters  was  outlined.  The  extreme  values  of  the  core  radius, 
in  that  region  are  3475.  7 and  3478.  8.  It  was  shown  that  the  solution 
should  be  sought  in  the  domain  of  the  transformed  parameters,  since  r 
and  dv/dz  are  closely  correlated. 


Onofrash,  N !.  , A.  V.  Drumya,  and  A.  A. 

Roman.  Qualitative  estimate  of  seismic  effect 
using  descriptive  data.  AN  Mold  SSR.  Izvestiya, 
no.  2,  1974,  80-00. 

Five  qualitative  methods  are  proposed  for  mapping  isoseismal 
curves.  The  problem  was  formalized  by  methods  of:  1)  information  theory; 

2)  hypothesis  v erification;  3)  parameter  evaluation;  4)  image  recognition; 
and  3)  taxonomy.  Ihe  first  four  methods  are  based  on  use  of  all  information 
contained  in  factual  material  as  well  as  prior  information  and  intensity  scale. 
The  last  method  is  based  exclusively  on  factual  material.  J he  methods  are 
illustrated  by  examples. 


1 


Bardan,  V.  Seismic  trace  resampling.  Studii 

si  cercetari  de  geologic,  geofizica,  geografie. 

Geofizica,  v.  12,  1974,  97-l',6. 

In  the  digital  recording  of  seismic  data,  the  sampling 
interval  is  1,  2,  4 or  8 msec,  depending  upon  the  highest  frequency  of  the 
seismic  signal.  Often,  to  increase  the  speed  of  data  processing,  a changing 
of  the  sampling  rate  is  used  (for  instance  from  2 msec  to  4 msec).  This 
paper  is  an  attempt  to  present  the  resampling  theory  of  seismic  data. 

In  the  first  part  of  the  article  sampling  theory  and  aliasing 
phenomenon  are  presented.  The  method  for  resampling  the  seismic  trace 
with  a different  interval  is  the  content  of  the  second  part.  The  design  of  a 
digital  minimum  phase  antialiasing  filter  is  presented  in  the  last  part. 


Iosif,  T. , and  S.  Iosif.  Optimization  of  seismic 
station  distribution  in  Rumania.  Studii  si  cercetari 
de  geologie,  geofizica,  geografie.  Geofizica,  v.  12, 

1974,  51-88. 

In  order  to  analyze  the  distribution  of  the  seismic  stations  in 
Romania,  standard  errors  of  the  earthquakes  parameters  (origin  time  ard 
focus  position)  for  that  territory  have  been  calculated.  The  values  of  the 
standard  errors  depend  on  the  spatial  distribution  of  the  seismic  stations 
ar.d  can,  therefore,  be  used  to  make  quantitative  comparison  among  different 
distributions  of  the  stations. 

Nine  combinations  of  stations,  including  the  existing  network 
and  proposed  locations  for  new  seismic  stations,  are  described  (see  Fig.  1 and 
Table  1).  If  the  locations  of  new  seismic  stations  are  properly  chosen,  the 
accuracy  in  determing  focal  parameters  improves  significantly.  The 
calculations  used  for  the  study  are  based  on  the  Monte  Carlo  method. 


-64- 


nucurejli  IT21.S  26*05.8' 

Clmpiilung  45  10,1  25*02,3 
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7.  Chela  ; 45  28  25*57' 

| 45*53'  22*55 

8 Timisoara  15*13,0'  21*13,5 


Galit/in 
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46*45' 

23*39' 

13. 

41*12 

28*38' 

14. 

47*40' 

23*30 

15. 

44*15 

23*15' 
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- NS.  EW 

VKGIK 

- NS,  EW,  Z 
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- NS.  EW,  Z 

Mccanic 
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VEGIK 

NS.  EW,  Z 

Klrnos 

7, 

Mccanic 

150  kg 

- NS,  EW 

VEGIK 

- NS,  EW, 
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VEGIK 

- NS.  EW, 

Z 

Klrnos 

- NS.  EW, 
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Mccanic 

450  kg 

- NS.  EW 

200  kg 

- NS.  EW 

VEGIK 

- NS,  EW, 

Z 

SKM 

- NS,  EW, 

Z 

Table  1.  Rumanian  Seismic  Stations. 


Malischewsky , P.  1 he  influence  of  curved 
discontinuities  on  the  propagation  of  seismic 
surface  waves  (In  English).  Gerlands  beitrage 
7.ur  geophysik,  v.  83,  no.  "1,  1974,  355-362. 

The  interaction  of  seismic  surface  waves  with  curved  discontinuities 
is  considered.  This  is  realized  for  the  first  time  on  the  basis  of  Alsop’s 
method  by  using  a curvilinear  co-ordinate  system.  The  obtained  reflection 
and  transmission  coefficients  are  complex;  that  means,  the  corresponding 
phase  shifts  of  seismic  surface  waves  at  curved  discontinuities  can  be  different 

from  0 or 
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Bisztricsany,  E. , and  Gy.  Szeidovitz. 

Strain  seismograph  (In  English).  Acta 
geodaetica,  geophysica  et  montanistica, 
v.  8,  nos.  3-4,  1973,  483-488. 

At  the  Budapest  Seismological  Observatory  the  first 
Hungarian  strain-seismometer  has  been  built.  The  22  m long  strainmeter 
has  been  constructed  from  55  porcelain  rods  of  65  mm  diameter,  thus  its 

cost  is  rather  low.  The  thermal  expansion  coefficient  of  the  porcelain 
used  is  not  much  greater  than  that  of  the  quartz. 


Guterch,  A.  Refraction  studies  of  structure  of 
the  Earth's  crust  and  upper  mantle  with  deep 
seismic  sounding  method  on  the  territory  of  Poland 
(In  English).  Acta  geophysica  Polonica,  v.  22, 
no.  3,  1974,  225-246. 

The  paper  presents  consecutive  stages  and  major  results  of 
investigations  of  the  earth's  crust  and  upper  mantle  made  with  the  deep 
seismic  sounding  method  in  Poland  in  the  years  1961-1973.  Basic  properties  of 
the  deepcrustal  structure  have  been  determined  along  the  SW-NE  direction  inter- 
secting the  main  tectonic  units  of  Poland.  Numerous  deep  fractures  were 
found  to  occur;  they  are  of  great  importance  in  the  tectonics  of  Poland,  in 
particular  in  the  problem  of  determining  the  margin  of  the  pre-Cambrian 
platform  of  Eastern  Europe,  which  is  one  of  the  major  tectonic  problems 
in  Europe.  The  structure  and  physical  properties  of  the  transition  zone 
between  the  crust  and  upper  mantle  have  been  determined  in  the  fore-Sudetic 
region  and  the  zone  of  the  Teisseyre-Tornquist  line. 


Chen  Pei-shau  and  Yen  Shou-min.  Relation 
between  seismic  «ourcc  mechanism  and  intensity 
distribution  (In  Chinese).  Acta  geophysica  Sinica, 
v.  18,  no.  1,  1975,  11-25. 

In  this  paper,  the  authors  combine  the  Haskell's  model  which 
is  a moving  fault  of  finite  length  and  the  Aki's  so-called  UZ  model  which  is 
constructed  by  fitting  an  exponentially  decaying  function  to  the  autocorrelation 
function  of  the  dislocation  velocity.  Neglecting  some  secondary  factors 
and  terms  and  at  the  same  time  considering  its  relation  to  earthquake 
intensity,  the  authors  have  tjoposed  a model  for  computing  the  distribution 
earthquake  intensity.  Using  this  model  and  the  seismic  source  parameters 
obtained  by  source  mechanism  studies,  the  authors  computed  the  theorectical 
distributions  of  earthquake  intensity  for  four  typical  earthquakes.  Comparison 
of  the  theoretical  with  the  observed  distributions  shows  the  agreement  is  close. 
It  shows  that  this  model  not  only  can  explain  many  earthquake  intensity 
distributions  in  practice,  some  of  which  are  hard  to  explain  from  the  point  of 
view  of  geological  structures,  but  also  may  be  used  for  computing  the 
theoretical  earthquake  intensity  distribution  for  a very  large  area. 


Kao  Lung-shen,  and  Ge  Huan-chen.  A preliminary 
study  of  P-  and  S-wave  velocities  under  high  pressure 
for  rock  samples  from  the  mainland  of  China  (In 
Chinese).  Acta  geophysica  Sinica,  v.  18,  no.  1,  1975, 

26-38. 

The  compressional  (P)  and  transversal  (S)  wave  velocities  of 
20  specimens  of  10  different  rock  types  and  1 specimen  of  silver  chloride  were 
determined  under  high  pressure.  The  pressure  in  measuring  the  P-wave 
velocity  of  three  of  the  specimens  went  up  as  high  as  20,  000  kg/cm2.  It  was 
found  that  above  4,000  kg/cm2  of  pressure,  the  change  of  velocity  with  pressure 
became  practically  linear,  but  the  rate  of  change  decreased  gradually  when  the 
pressure  was  greater  than  10,  000  kg /cm2. 
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For  the  determination  of  the  S-wave  velocity,  in  order  to  obtain 
more  accurate  first  arrival  times,  coupled  circuits  were  connected  respectively 
to  both  the  emitting  and  receiving  crystals.  By  employing  a method  of 
"extrapolation  of  multiple  points"  measurement  of  the  S-wave  train,  we  can 
reasonably  raise  the  accuracy.  Errors  of  the  results  were  then  analysed. 


Liu  Yi-ming,  Yao  Hung,  and  Chou  Hai-ian. 

Analysis  of  the  characteristics  and  trend  of 
seismic  activities  in  the  Keping,  rupture  zone 
in  Sinkiang,  China  (In  Chinese).  Acta  geophysica 
Sinica,  v.  18,  no.  1,  1975,  39-51. 

Using  the  results  from  the  analysis  of  the  Keping  earthquake 
of  6.2  magnitude  (Richter  scale)  in  western  Sinkiang  in  16  January  1972,  a 
further  study  of  the  seismic  activities  for  a larger  region  around  Keping  is 
made.  The  characteristics  of  seismicity,  source  mechanisms,  velocity 
variations,  as  well  as  tectonic  movements  of  the  region  are  correlated  and 
analyzed.  Basing  upon  the  results  obtained,  a preliminary  discussion  of  the 
trend  of  seismic  activities  in  the  region  is  attempted. 
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